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The membrane paradigm is the remarkable view that, to an external observer, 



t^. ' a black hole appears to behave exactly like a dynamical fluid membrane, 

q ' obeying such pre-relativistic equations as Ohm's law and the Navier-Stokes 

^ ' equation. It has traditionally been derived by manipulating the equations 

of motion. Here we provide an action formulation of this picture, clarifying 
what underlies the paradigm, and simplifying the derivations. Within this 



framework, we derive previous membrane results, and extend them to dy- 
onic black hole solutions. We discuss how it is that an action can produce 
dissipative equations. Using a Euclidean path integral, we show that famil- 
iar semi-classical thermodynamic properties of black holes also emerge from 
the membrane action. Finally, in a Hamiltonian description, we establish the 
validity of a minimum entropy production principle for black holes. 
PACS: 03.50.-z, 04.20.-q, 04.70.-s, 04.70.Bw, 04.70.Dy 
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I. INTRODUCTION 



The event horizon of a black hole is a peculiar object: it is a mathematically defined, 
locally undetectable boundary, a surface-of-no-return inside which light cones tip over and 
"time" becomes spatial ffl]. Otherwise natural descriptions of physics often have trouble 
accommodating the horizon; as the most primitive example, the familiar Schwarzschild 
metric has a co-ordinate singularity there. Theories of fields that extend to the horizon 
face the additional challenge of having to define boundary conditions on a surface that is 
infinitely red-shifted, has a singular Jacobian, and possesses a normal vector which is also 
tangential. These considerations might induce one to believe that black hole horizons are 
fundamentally different from other physical entities. 

On the other hand, further work has established a great variety of analogies between the 
horizon and more familiar, pre-relativistic bodies. In addition to the famous four laws of 
black hole thermodynamics 0-§J, which are global statements, there is also a precise local 
mechanical and electrodynamic correspondence. In effect, it has been shown |6|-|T0|j that 
an observer who remains outside a black hole perceives the horizon to behave according 
to equations that describe a fluid bubble with electrical conductivity as well as shear and 
bulk viscosities. Moreover, it is possible to define a set of local surface densities, such as 
charge or energy-momentum, which inhabit the bubble surface and which obey conservation 
laws. Quite remarkably, a general-relativistically exact calculation then leads, for arbitrary 
non-equilibrium black holes, to equations for the horizon which can be precisely identified 
with Ohm's law, the Joule heating law, and the Navier-Stokes equation. 

These relations were originally derived for the mathematical, or true, event horizon. For 
astrophysical applications it became more convenient to consider instead a "stretched hori- 
zon," a 2+1- dimensional time-like surface located slightly outside the true horizon. Because 
it has a non-singular induced metric, the stretched horizon provides a more tractable bound- 
ary on which to anchor external fields; outside a complicated boundary layer, the equations 
governing the stretched horizon are to excellent approximation [PIPI the same as those for 
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the true horizon. This view of a black hole as a dynamical time-like surface, or membrane, 



has been called the membrane paradigm [13 



Most of the mentioned results have been derived through general-relativistic calculations 
based on various intuitive physical arguments. In this paper, we show that the gravitational 
and electromagnetic descriptions of the membrane can be derived systematically, directly, 
and more simply from the Einstein-Hilbert or Maxwell actions. Aside from the appeal 
inherent in a least action principle, an action formulation is a unifying framework which is 
easily generalizable and has the advantage of providing a bridge to thermodynamics and 



quantum mechanics (see |14| for related work). In a follow-up paper, we exploit these 
advantages to evaluate some effects of back-reaction of spacetime geometry on Hawking 
radiation fl5fl . 

The key idea in what follows is that, since (classically) nothing can emerge from a black 
hole, an observer who remains outside a black hole cannot be affected by the dynamics inside 
the hole. Hence the equations of motion ought to follow from varying an action restricted to 
the external universe. However, the boundary term in the derivation of the Euler-Lagrange 
equations does not in general vanish on the stretched horizon as it does at the boundary of 
spacetime. In order to obtain the correct equations of motion, we must add to the external 
action a surface term that cancels this residual boundary term. The membrane picture 
emerges in interpreting the added surface term as electromagnetic and gravitational sources 
residing on the stretched horizon. 

In the rest of this paper, we examine individually the boundary terms for the electro- 
magnetic, gravitational, and axidilaton cases. We also discuss dissipation and provide ther- 
modynamic and Hamiltonian descriptions. We use lowercase indices for four- dimensional 
tensor indices and uppercase indices for the two-dimensional tensors that occupy space-like 
sections of the horizon. We use geometrized units (G = c = 1), and a spacetime metric 



with positive signature ( — h ++). Our sign conventions are those of MTW |T6[, with the 
exception of the extrinsic curvature which we define to have positive trace for a convex 
surface. 



II. HORIZON PRELIMINARIES 



In this section, we fix our conventions, first in words, then in equations. Through every 
point on the true horizon there exists a unique null generator l a which we may parameterize 
by some regular time co-ordinate whose normalization we fix to equal that of time-at- infinity. 
Next, we choose a time-like surface just outside the true horizon. This is the stretched 
horizon, 7i, whose location we parameterize by a <C 1 so that a — » is the limit in which 
the stretched horizon coincides with the true horizon. We will always take this limit at the 
end of any computation. Since many of the useful intermediate quantities will diverge as 
inverse powers of a, we renormalize them by the appropriate power of a. In that sense, a 
plays the role of a regulator. 

For our purposes, the principal reason for preferring the stretched horizon over the true 
horizon is that the metric on a time-like - rather than null - surface is non-degenerate, 
permitting one to write down a conventional action. Generically (in the absence of horizon 
caustics), a one-to-one correspondence between points on the true and stretched horizons 



is always possible via, for example, ingoing null rays that pierce both surfaces (see |L2| for 
details). 

We can take the stretched horizon to be the world-tube of a family of time-like observers 
who hover just outside the true horizon. These nearly light-like "fiducial" observers are 
pathological in that they suffer an enormous proper acceleration and measure quantities 
that diverge as a — > 0. However, although we take the mathematical limit in which the true 
and stretched horizons conflate, for physical purposes the proper distance of the stretched 
horizon from the true horizon need only be smaller than the length scale involved in a 
given measurement. In that respect, the stretched horizon, although a surrogate for the 
true horizon, is actually more fundamental than the true horizon, since measurements at 
the stretched horizon constitute real measurements that an external observer could make 
and report, whereas accessing any quantity measured at the true horizon would entail the 
observer's inability to report back his or her results. 
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We take our fiducial observers to have world lines U a , parameterized by their proper time, 
r. The stretched horizon also possesses a space-like unit normal n a which for consistency we 
shall always take to be outward-pointing. Moreover, we choose the normal vector congruence 
on the stretched horizon to emanate outwards along geodesies. We define a by requiring 
that aU a — > l a and an a — > l a ; hence aU a and an" are equal in the true horizon limit. This 
is nothing more than the statement that the null generator l a is both normal and tangential 
to the true horizon, which is the defining property of null surfaces. Ultimately though, it 
will be this property that will be responsible for the dissipative behavior of the horizons. 
The 3-metric, h ab , on H can be written as a 4-dimensional tensor in terms of the spacetime 
metric and the normal vector, so that hi projects from the spacetime tangent space to the 
3-tangent space. Similarly, we can define the 2-metric, Jab, of the spacelike section of TC 
to which U a is normal, in terms of the stretched horizon 3-metric and U a , thus making 
a 2+1+1 split of spacetime. We denote the 4-covariant derivative by V a , the 3-covariant 
derivative by \ a , and the 2-covariant derivative by p. For a vector in the stretched horizon, 
the covariant derivatives are related by h c d V c w a = w" d — K c d w c n a where K b = h^V ' c n a is the 
stretched horizon's extrinsic curvature, or second fundamental form. In summary, 



(2.1) 




U 2 = -1 , lim aU a = I 



a 



(2.2) 



n 2 = +1 , a c = n a V a n c = , lim an a = l a 



(2.3) 



h a h = g a b - n a n b , ^ = h a b + U a U b = g a b - n a n b + U a U b 



(2.4) 



K a b =h c b W c n\ K ab = K ba , K ab n b = 



(2.5) 



w c eH => h c d V c w a = 



W\ d - K d w c n a =^ V c w c = W\ c + w c a c = 




(2.6) 



5 



The last expression relates the covariant divergence associated with g ab to the covariant 
divergence associated with h ab . 

For example, the Reissner-Nordstrom solution has 

ds 2 = - il - — + dt 2 + ll- — + -^\ dr 2 + r 2 dfl 2 , (2.7) 

so that a stretched horizon at constant r would have 
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a =( 1 " — + ' ^ 

?7 a = -a(dt) a , (2.9) 

and 

n a = +a' 1 (dr) a . (2.10) 



III. THE ACTIONS 

To find the complete equations of motion by extremizing an action, it is not sufficient to 
set the bulk variation of the action to zero: one also needs to use the boundary conditions. 
Here we take our Dirichlet boundary conditions to be 5(f) = at the singularity and at the 
boundary of spacetime, where (j) stands for any field. 

Now since the fields inside a black hole cannot have any classical relevance for an external 
observer, the physics must follow from varying the part of the action restricted to the 
spacetime outside the black hole. However, this external action is not stationary on its own, 
because boundary conditions are fixed only at the singularity and at infinity, but not at the 
stretched horizon. Consequently, we re-write the total action as 

Sworid = \Sout + S sur f) + (Si n — S sur f) , (3-1) 

where now SS out + SS sur f = 0, which implies also that SSi n — SS sur f = 0. The total action 
has been broken down into two parts, both of which are stationary on their own, and which 
do not require any new boundary conditions. 
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The surface term, S sur f, corresponds to sources, such as surface electric charges and 
currents for the Maxwell action, or surface stress tensors for the Einstein-Hilbert action. 
The sources are fictitious: an observer who falls through the stretched horizon will not 
find any surface sources and, in fact, will not find any stretched horizon. Furthermore, the 
field configurations inside the black hole will be measured by this observer to be entirely 
different from those posited by the membrane paradigm. On the other hand, for an external 
fiducial observer the source terms are a very useful artifice; their presence is consistent with 
all external fields. This situation is directly analogous to the method of image charges in 
electrostatics, in which a fictitious charge distribution is added to the system to implement, 
say, conducting boundary conditions. By virtue of the uniqueness of solutions to Poisson's 
equation with conducting boundary conditions, the electric potential on one - and only one 
- side of the boundary is guaranteed to be the correct potential. An observer who remains 
on that side of the boundary has no way of telling through the fields alone whether they 
arise through the fictitious image charges or through actual surface charges. The illusion is 
exposed only to the observer who crosses the boundary to find that not only are there no 
charges, but the potential on the other side of the boundary is quite different from what it 
would have been had the image charges been real. 



In the rest of this section, we shall implement eqn. |3J] concretely in important special 
cases. 

A. The Electromagnetic Membrane 

The external Maxwell action is 

Sout[A a ] = J d'x^g (-y^ 2 + J ■ A) , (3.2) 

where F is the electromagnetic field strength. Under variation, we obtain the inhomogeneous 
Maxwell equations, 

V b F ab = 4vrJ a , (3.3) 



as well as the boundary term 

— / d 3 xV^hF ab n a 5A b , (3.4) 
An J 

where h is the determinant of the induced metric, and n a is the outward-pointing space-like 
unit normal to the stretched horizon. We need to cancel this term. Adding the surface term 

S S urf[A a ] = +J d 3 xV^hj s ■ A , (3.5) 

we see that we must have 

1 



f s = +^F ah n b . (3.6) 

The surface 4-current, j®, has a simple physical interpretation. We see that its time- 
component is a surface charge, a, that terminates the normal component of the electric 
field just outside the membrane, while the spatial components, j s , form a surface current 
that terminates the tangential component of the external magnetic field: 

E± = -U a F ab n b = Ana (3.7) 

B( = ei^F ab n b = An (j s x h) A . (3.8) 

It is characteristic of the membrane paradigm that a and j s are local densities, so that the 
total charge on the black hole is the surface integral of a over the membrane, taken at some 
constant universal time. This is in contrast to the total charge of a Reissner-Nordstrom 
black hole, which is a global characteristic that can be defined by an integral at spatial 
infinity. 



From Maxwell's equations and eqn. 3J;, we obtain a continuity equation for the mem- 
brane 4-current which, for a stationary hole, takes the form 

|^ + V 2 -J S = -J\ (3.9) 

where V2 • j s = {la3s) A ls ^ ne two-dimensional divergence of the membrane surface current, 
and — J n = —J a n a is the amount of charge that falls into the hole per unit area per unit 
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proper time, r. Physically, this equation expresses local charge conservation in that any 
charge that falls into the black hole can be regarded as remaining on the membrane: the 
membrane is impermeable to charge. 

The equations we have so far are sufficient to determine the fields outside the horizon, 
given initial conditions outside the horizon. A plausible requirement for initial conditions at 
the horizon is that the fields measured by freely-falling observers (FFO's) at the stretched 
horizon be finite. There being no curvature singularity at the horizon, inertial observers who 
fall through the horizon should detect nothing out of the ordinary. In contrast, the fiducial 
observers (FIDO's) who make measurements at the membrane are infinitely accelerated. 
Their measurements, subject to infinite Lorentz boosts, are singular. For the electromagnetic 
fields we have, with 7 the Lorentz boost and using orthonormal co-ordinates, 

E! IDO « 7 (e[ FO - B F /°) , B$ ID ° « 7 {b ffo - Ef F °) , (3.10) 
B! IDO « 7 K F ° - E F /°) , Ef°° w 7 (eFFO _ b ffo^ (3 n) 

or, more compactly, 

E FIDO = nxB FIDO . (3.12) 
That is, the regularity condition states that all radiation in the normal direction is ingoing; 



a black hole acts as a perfect absorber. Combining the regularity condition with eqn. p7c 
and dropping the FIDO label, we arrive at 

E\\ = Anj s ■ (3.13) 

That is, black holes obeys Ohm's law with a surface resistivity of p = An rs 377 Q. Further- 
more, the Poynting flux is 

S = ±-(ExB) = -f sP h. (3.14) 

We can integrate this over the black hole horizon at some fixed time. However, for a generic 
stretched horizon, we cannot time-slice using fiducial time as different fiducial observers have 
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clocks that do not necessarily remain synchronized. Consequently we must use some other 
time for slicing purposes, such as the time at infinity, and then include in the integrand a 
(potentially position-dependent) factor to convert the locally measured energy flux to one at 
infinity. With a clever choice of the stretched horizon, however, it is possible to arrange that 
all fiducial observers have synchronized clocks. In this case, two powers of a, which is now 
the lapse, are included in the integrand. Then, for some given universal time, t, the power 
radiated into the black hole, which is also the rate of increase of the black hole's irreducible 
mass, is given by 

^- = -J^S-dA=+ ja 2 j 2 s pdA. (3.15) 

That is, black holes obey the Joule heating law, the same law that also describes the dissi- 
pation of an ordinary ohmic resistor. 



B. The Gravitational Membrane 

We turn now to gravity. The external Einstein-Hilbert action is 

Sout[g ab ] = -^~f d 4 x^gR+^-fd 3 xV±hK + S matter , (3.16) 

where R is the Ricci scalar, K is the trace of the extrinsic curvature, and where for conve- 
nience we have chosen the field variable to be the inverse metric g ab . The surface integral 
of K is only over the outer boundary of spacetime, and not over the stretched horizon. It 
is required in order to obtain the Einstein equations because the Ricci scalar contains sec- 
ond order derivatives of g a b- When this action is varied, the bulk terms give the Einstein 
equations, 

R ab -\gabR^^T ab . (3.17) 

We are interested however in the interior boundary term. This comes from the variation of 
the Ricci tensor. We note that 
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g ab 8R ab = V a ( V 6 (Sg ab ) - g cd V a (5g cd )) , (3.18) 

where 5g a b = —g ac gbd5g cd ■ Gauss' theorem now gives 

J d A x^~g (g ab 5R ab ) =-j d?x^hn a h bc (V c (5g ab ) - V a (5g bc )) , (3.19) 

where the minus sign arises from choosing n a to be outward-pointing. Applying the Leibniz 
rule, we can rewrite this as 

J d 4 x^g~ (g ab 5R ab ) = J d 3 xV^hh bc (V a (n a 5g bc ) -5g bc V a (n a ) 

- V c (n a 5g ab ) + 5g ab V c (n a )) . (3.20) 

Now, in the limit that the stretched horizon approaches the null horizon, the first and 
third terms on the right-hand side vanish: 

J d 3 xV^hh bc (V„ (n a 5g bc ) - V c (n a 5g ab )) = . (3.21) 

A proof of this identity is given in the appendix. With K ba = h bc 'S7 c n a ) the variation of the 
external action is 

SS out [g ab ] = -L [ d 3 xV^h (Kh ab - K ab ) Sg ab . (3.22) 

Since the expression in parentheses contains only stretched horizon tensors, the normal 
vectors in the variation Sg ab = 5h ab + 5n a n b + n a Sn b contribute nothing. As in the electro- 
magnetic case, we add a surface source term to the action to cancel this residual boundary 
term. The variation of the required term can therefore be written as 

5S mrf [h ab ] = ~ / d 3 xV^ht sab 5h ab . (3.23) 

We shall see later that this variation is integrable i.e. an action with this variation exists. 
Comparison with equation |3.22j yields the membrane stress tensor: 



tf = +— (Kh ab -K ab ) . (3.24) 
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Now just as a surface charge produces a discontinuity in the normal component of the electric 
field, a surface stress term creates a discontinuity in the extrinsic curvature. The relation 



between the discontinuity and the source term is given by the Israel junction condition 16 



tf = ^([K]h ab -[Kr) , (3.25) 

where [K] = K + — K- is the difference in the extrinsic curvature of the stretched horizon 
between its embedding in the external universe and its embedding in the spacetime internal 



to the black hole. Comparing this with our result for the membrane stress tensor, eqn. |3.24j , 
we see that 

K ab = , (3.26) 
so that the interior of the stretched horizon molds itself into flat space. The Einstein 



equations, eqn. 3.17, can be rewritten via the contracted Gauss-Codazzi equations [16 



as 



tf ]b = -h a c T cd n d . (3.27) 



Eqns. |3.24 and |3.27| taken together imply that the stretched horizon can be thought of 



as a fluid membrane, obeying the Navier-Stokes equation. To see this, recall that as we send 
a to zero, both aU a and an a approach l a , the null generator at the corresponding point on 
the true horizon. Hence, in this limit we can equate aU a and an a , permitting us to write the 
relevant components of Kg, in terms of the surface gravity, g, and the extrinsic curvature, 
kg, of a spacelike 2-section of the stretched horizon: 

U c V c n a -> a- 2 l c V c l a = a- 2 g H l a => K b a U a U b = -g = -a^gn, K§ = it K bU b = , (3.28) 
where gn = otg is the renormalized surface gravity at the horizon, and 

7.Vy - v/ K» = l a A K h al b B = a~ l k B A , (3.29) 

where )zab is the extrinsic curvature of a spacelike 2-section of the true horizon, 

kAB = lA l B\\d = ^£l*lAB , (3.30) 
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where £i» is the Lie derivative in the direction of l a . We can decompose k^B into a traceless 
part and a trace, kAB = &ab + \iab9, where gab is the shear and the expansion of the 
world lines of nearby horizon surface elements. Then 

tf = -L f-a AB + 7 AB ( U + g]) . (3.31) 



8tt V V2 

But this is just the equation for the stress of a two-dimensional viscous Newtonian fluid 
fT7fl with pressure p = f-, shear viscosity rj = ■A-, and bulk viscosity ( = — ■A-. Hence we 
may identify the horizon with a two-dimensional dynamical fluid, or membrane. Note that, 
unlike ordinary fluids, the membrane has negative bulk viscosity. This would ordinarily 
indicate an instability against generic perturbations triggering expansion or contraction. It 
can be regarded as reflecting a null hypersurface's natural tendency to expand or contract 
lOfl . Below we shall show how for the horizon this particular instability is replaced with a 



different kind of instability. 

Inserting the A-momentum density tJ^j^Ub = t^ A = tta into the Einstein equations, 



eqns. |3.27| , we arrive at the Navier-Stokes equation, 

£tKa = -VaP + C^aO + 2 W B {{B - T2 , (3.32) 

where £ t iia — jff is the Lie derivative (which is the general-relativistic equivalent of the 
convective derivative) with respect to proper time, and — T\ = —^i a A T^n c is the flux of 
A-momentum into the black hole. 

Inserting the [/-momentum (energy) density t^ b U a U b = S = — ^ gives 

£ r S + OH = -p6 + CO 2 + 2 WA b(t ab + T b a n a U b , (3.33) 



which is the focusing equation for a null geodesic congruence [18]. We might now suspect 



that if the analogy with fluids extends to thermodynamics then eqn. 3.33, as the equation 



of energy conservation, must be the heat transfer equation |T7] for a two-dimensional fluid. 



Writing the expansion of the fluid in terms of the area, AA, of a patch, 



1 dAA 
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we see that we can indeed re-write eqn. [3.33| as the heat transfer equation (albeit with an 
additional relativistic term on the left), 

T (dAS_ = f CQ 2 + 2t)GabG ab + T a naU ^ AA j (3>35) 



\ dr g dr 2 J 

with T the temperature and S the entropy, provided that the entropy is given by 



and the temperature by 



S = r&A , (3.36) 



T = ^4—9 , (3.37) 



where r/ is some proportionality constant. 

Thus, the identification of the horizon with a fluid membrane can be extended to the 
thermodynamic domain. Nonetheless, the membrane is an unusual fluid. The focusing equa- 
tion itself, eqn. |3.33j , is identical in form to the equation of energy conservation for a fluid. 



However, because the energy density, E, is proportional to the expansion, S = — one ob- 
tains a non-linear first-order differential equation for 9 which has no counterpart for ordinary 
fluids. The crucial point is that, owing to the black hole's gravitational self-attraction, the 
energy density is negative, and the solution to the differential equation represents a horizon 
that grows with time. For example, the source-free solution with a time-slicing for which 
the horizon has constant surface gravity is 



Because of the sign of the exponent, this would represent an ever-expanding horizon if 9 (to) 
were an initial condition; the area of the horizon, which is related to 6 by 9 = j- In ^ , 
expands exponentially with time. To avoid this runaway, one must impose "teleological 
boundary conditons" (that is, final conditions) rather than initial conditions. Hence, the 
horizon's growth is actually acausal; the membrane expands to intercept infalling matter 
that is yet to fall in ||13|| . This is because the membrane inherits the global character of 
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the true horizon: the stretched horizon covers the true horizon whose location can only 
be determined by tracking null rays into the infinite future. In fact, the left-hand side of 
the heat transfer equation, eqn. |3.35| , is of the same form as that of an electron subject to 
radiation reaction; the acausality of the horizon is therefore analogous to the pre-acceleration 
of the electron. 

At this classical level, using only the equations of motion, the parameter r] in eqn. 



3.36| is undetermined. However, since we have an action, we hope to do better, since the 
normalization in the path integral is now fixed. By means of a Euclidean path integral, we 
should actually be able to derive the Bekenstein-Hawking entropy, including the coefficient 
rj, from the membrane action. We do this in a later section. 

C. The Axidilaton Membrane 

Another advantage of the action formulation is that it is easily generalized to arbitrary 
fields. For example, we can extend the membrane paradigm to include the basic fields of 
quantum gravity. Here we use the tree-level effective action obtained from string theory 
after compactification to four macroscopic dimensions. This action is a generalization of 
the classical Einstein-Hilbert-Maxwell action to which it reduces when the axidilaton, A 
(sometimes written as S), is set to jj-. The action is 

S[X,X,A a ,g ab ] = J d'x^g (JL - |^ + i (Ai? - XF 2 _ )) , (3.39) 

where R is the four-dimensional Ricci curvature scalar, F± = F ± iF are the self- and anti- 
self-dual electromagnetic field strengths, and A = Ai +i\>z = a + ie~ 2< ^ is the axidilaton, with 
a the axion and <p the dilaton. Solutions to the equations of motion arising from this action 



include electrically- (Reissner-Nordstrom) and magnetically-charged black holes [115,20], as 
well as their duality-rotated cousins, dyonic black holes pl| , which carry both electric and 
magnetic charge. 

The equations of motion are 
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„ fd a X\ \d\\ 2 i 2 , in . 

and 

V a (AFf - XF^ b ) = , (3.41) 

besides the Einstein equations. 

As before, we require the external action to vanish on its own. Integration by parts on 
the axidilaton kinetic term leads to a variation at the boundary, 

where n a is again chosen to be outward-pointing. To cancel this, we add the surface term 

S surf = J d 3 xV^h (\q + Ag) , (3.43) 

so that 



q = — . (3.44) 

A 2 



To interpret this, we note that the kinetic term in A is invariant under global SL (2, M) 
transformations of the form 

x ^aX + b ad _ bc=1 (3 >45) 
cX + d 

which are generated by Peccei-Quinn shifts, Ai — > Ai + b, and duality transformations, 
A — > — The Peccei-Quinn shift of the axion can be promoted to a classical local symmetry 
to yield a Nother current: 

Jp- Q = -^(d a X + d a X) . (3.46) 
Therefore, under a Peccei-Quinn shift, 



f = J d 3 xV^h5X(q + q) = J d 3 xV^h5X(n a J a P _ Q ) . (3.47) 
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The sum of the q and q terms induced at the membrane, eqn. |3.44| , is the normal component 
of the Peccei-Quinn current. Hence, at the membrane, 



(kj: 



b 

P-Q 



| a 



•FF — V a ((q + q) n a ) 



(3.48) 



That is, the membrane term V a {{q + q) n a ) augments the dyonic FF term as a source for 
the three-dimensional Peccei-Quinn current, h^J P _Q, at the membrane. 

The membrane is again dissipative with the Peccei-Quinn charge accounting for the 
dissipation in the usual a->0 limit. The local rate of dissipation is given by the bulk stress 
tensor at the membrane: 



2|„|2 



rp TT a^b_ 1 WA + WA t \ 2 2 \q 



(3.49) 
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In addition, the presence of the axidilaton affects the electromagnetic membrane. (The 
gravitational membrane is unaffected since the surface terms come from the Ricci scalar 
which has no axidilaton factor.) The electromagnetic current is now 



?ab XF^ b ) n b 



f s = -2i (XFl 
The surface charge is therefore 

a = 4 (X 2 E ± + X 1 B ± ) , 

and the surface current is 

j s = 4 (x 2 n x B\\ — X\h x E\^j , 



(3.50) 



(3.51) 



(3.52) 



which, by the regularity of the electromagnetic field, eqn. |3.12|, satisfies 
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(3.53) 



The conductivity is now a tensor. When the axion is absent, the resistivity is 

1 



P 



4A 2 ' 
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(3.54) 



The inverse dependence on A2 is to be expected on dimensional grounds. The pure dilaton 
action can be derived from Kaluza-Klein compactification of pure gravity in five dimensions, 
where the fifth dimension is curled into a circle of radius e~ 2 ^ = A2. In five dimensions, 
with c = 1, resistance (and hence resistivity for a two-dimensional resistor such as the 
membrane) has dimensions of inverse length. Using the regularity condition, eqn. |3.12| , the 
rate of dissipation, for a stretched horizon defined to have uniform lapse a with respect to 
time at infinity, t, is 

d J^L = _J a 2s.dA = J Aa 2 \ 2 E 2 dA = J a 2 ^-J s 2 dA , (3.55) 

which is the Joule heating law in the presence of an axidilaton. 



IV. DISSIPATION 

Given that the bulk equations of motion are manifestly symmetric under time- reversal, 
the appearance of dissipation, as in Joule heating and fluid viscosity, might seem mysterious, 
all the more so since it has been derived from an action. 

The procedure, described here, of restricting the action to some region and adding surface 
terms on the boundary of the region cannot be applied with impunity to any arbitrary 
region: a black hole is special. This is because the region outside the black hole contains its 
own causal past; an observer who remains outside the black hole is justified in neglecting 
(indeed, is unaware of) events inside. However, even "past sufficiency" does not adequately 
capture the requirements for our membrane approach. For instance, the past light cone of 
a spacetime point obviously contains its own past, but an observer in this light cone must 
eventually leave it. Rather, we define the notion of a future dynamically closed set: 

A set S in a time-orientable globally hyperbolic spacetime (M, g a b) is future dy- 
namically closed if J~{S) = S, and if, for some foliation of Cauchy surfaces 
£t parameterized by the values of some global time function, we have that 
Vt Vpe (SHEJ \/(t>t )3qe ( 1+ (p) n S n E t ) . 
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That is, S is future dynamically closed if it contains its own causal past and if from every 
point in S it is possible for an observer to remain in S. Classically, the region outside the 
true horizon of a black hole is dynamically closed. So too is the region on one side of a null 
plane in flat space; this is just the infinite-mass limit of a black hole. The region outside 
the stretched horizon is strictly speaking not dynamically closed since a signal originating 
in the thin region between the stretched horizon and the true horizon can propagate out 
beyond the stretched horizon. However, in the limit that the stretched horizon goes to the 
true horizon, a — > 0, this region becomes vanishingly thin so that in this limit, which is in 
any case assumed throughout, we are justified in restricting the action. 

The breaking of time-reversal symmetry comes from the definition of the stretched hori- 
zon; the region exterior to the black hole does not remain future dynamically closed un- 
der time-reversal. In other words, we have divided spacetime into two regions whose dy- 
namics are derived from two different simultaneously vanishing actions, 5 (S out + S sur f) = 
^ (Si n — Ssurf) = 0. Given data on some suitable achronal subset we can, for the exterior 
region, predict the future but not the entire past, while, inside the black hole, we can "post- 
diet" the past but cannot determine the entire future. Thus, our choice of the horizon as 
a boundary implicitly contains the irreducible logical requirement for dissipation, that is, 
asymmetry between past and future. 

Besides the global properties that logically permit one to write down a time-reversal 
asymmetric action, there is also a local property of the horizon which is the proximate cause 
for dissipation, namely that the normal to the horizon is also tangential to the horizon. 
Without this crucial property - which manifests itself as the regularity condition, or the 
identification of the stretched horizon extrinsic curvature with intrinsic properties of the 
true horizon - there would still be surface terms induced at the stretched horizon, but no 
dissipation. 

The regularity condition imposed at the boundary is not an operator identity, but a 
statement about physical states: all tangential electromagnetic fields as measured by a 
fiducial observer must be ingoing. Such a statement is not rigorously true. For any given 
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value of a — (l — ^) 2 , there is a maximum wavelength, X max , for outgoing modes that are 
invisible to the observer: 

r-2M , . 

Xmax = r -> 2Ma . (4.1) 

Dissipation occurs in the membrane paradigm because the finite but very high-frequency 
modes that are invisible to the fiducial observer are tacitly assumed not to exist. The reg- 
ularity condition amounts to a coarse-graining over these modes. It is conceivable that for 
a theory with benign ultraviolet behavior, the amount of information lost is finite. Einstein 
gravity is not such a theory, but one may ask abstractly whether an effective horizon theory 



could exist at a quantum level p2|j23[| . Quantum effects cause the black hole to emit radia- 
tion. In order to preserve time-evolution unitarity, we might require the emitted radiation to 
be correlated with the interior state of the black hole. In this case, the membrane viewpoint 
remains valid only as a classical description, since quantum-mechanically the external uni- 
verse receives information from the black hole in the form of deviations of the radiation from 
thermality; the crucial premise that the outside universe is emancipated from the internal 
state of the black hole is violated. It is important to emphasize, however, that correlations 
between the radiation and the horizon itself do not preclude the membrane paradigm. In- 
deed, the fact that the Bekenstein-Hawking entropy is proportional to the surface area of 
the black hole suggests that, even at the quantum level, an effective horizon theory may not 
be unfeasible. 

V. THE THERMODYNAMIC MEMBRANE 

To make contact with thermodynamics, we perform an analytic continuation to imagi- 
nary time, t = it, so that the path integral of the Euclideanized action becomes a partition 
function. For a stationary hole, regularity (or the removal of a conical singularity) dictates 



a period (3 = J dr = in imaginary time [23], where gn is the surface gravity; for a 



Schwarzschild hole, (3 = 8nM. This is the inverse Hawking temperature in units where 
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h = c = G = kB = l- The partition function is then the path integral over all Euclidean 
metrics which are periodic with period ^ in imaginary time. Since the dominant contri- 
bution to the path integral comes from the classical solution, we can evaluate the partition 
function in a stationary phase approximation: 

Z = J Dgfexp (-1 (S* ut {g% b } + Sf urf [h ab })) « exp (-1 (^[g&J + Sf urf [h a E b d })) • (5-1) 

The external action itself can be written as S out = Sbuik + <Soo> where Sbuik is zero for 
a black hole alone in the universe. The boundary term Soo is the integral of the extrinsic 
curvature of the boundary of spacetime. In fact, a term proportional to the surface area 
at infinity can be included in Soo without affecting the Einstein equations since the metric 
is held fixed at infinity during variation. In particular, the proportionality constant can be 
chosen so that the action for all of spacetime is zero for Minkowski space: 

S oo = ^~ [ d 3 xV^h[K] , (5.2) 

where [K] is the difference in the trace of the extrinsic curvature at the spacetime boundary 
for the metric g a b and the flat-space metric r] a b. With this choice, the path integral has 
a properly normalized probabilistic interpretation. The Euclideanized value of for the 
Schwarzschild solution is then |24] 



Sg = hm i- (-32tt 2 m) I (2r - 3M) -2r(l- ! I = +4ttM 2 . (5.3) 



To obtain an explicit action for the membrane, we must integrate its variation, eqn. 3.23: 



SS surf [h ab ] = - JL f d 3 xv^h (Kh ab - K ab ) 5h ab . (5.4) 

1D7T J 

We see that 

S S urf[h ab ] = J d 3 xV^h (B ab h ab - b) , (5.5) 

is a solution, provided that the (undifferentiated) source terms are B ab = +j^K ab and 
b = —r^-K. This action has the form of surface matter plus a negative cosmological constant 
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in three dimensions. The value of the membrane action for a solution to the classical field 
equations is then 

SsurAKi] = +F" / d 3 xyTh^K cl . (5.6) 
on J v 

To evaluate this, we can take our fiducial world-lines U a to be normal to the isometric 
time-slices of constant Schwarzschild time. The stretched horizon is then a surface of con- 
stant Schwarzschild r. Hence a = (l — ^jr) 2 , — 0, and K = g + 6 = g, the unrenormalized 



surface gravity of the stretched horizon. Inserting these into eqn. |5.6| , we obtain for the Eu- 
clidean action 

Sf urf = lim — ( ( -dr) aAixr 2 g = -nr H 2 = -4vrM 2 , (5.7) 

where r# = 2M is the black hole's radius, and gu = otg = -A* is its renormalized surface 
gravity. 

The Euclidean membrane action exactly cancels the external action, eqn. f>.3[ Hence the 
entropy is zero! That, however, is precisely what makes the membrane paradigm attractive: 
to an external observer, there is no black hole - only a membrane - and so neither a gen- 
eralized entropy nor a strictly obeyed second law of thermodynamics. The entropy of the 
outside is simply the logarithm of the number of quantum states of the matter outside the 
membrane. This number decreases as matter leaves the external system to fall through and 
be dissipated by the membrane. When all matter has fallen into the membrane, the outside 
is in a single state - vacuum - and has zero entropy, as above. 

To recover the Bekenstein-Hawking entropy, we must then use not the combination of 
external and membrane actions, which gave the entropy of the external system, but the 
combination of the internal and membrane actions, 

Z B -h = J Dgfexp (sf n [gf) - S? urf [h<$})) , (5.8) 

where now S sur f is subtracted (see eqn. |3.1p . With S in = J d^Xy/^g R = 0, the partition 
function for a Schwarzschild hole in the stationary phase approximation is 
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Z B _ H ^exp(-^(+47rM 2 )) , (5.9) 

from which the Bekenstein-Hawking entropy, Sb-h, immediately follows: 

Sb-h = p[m+ 1 ^j z A = SnM [M - ^4ttM 2 ) = , (5.10) 

which is the celebrated result. 

For more general stationary (Kerr-Newman) holes, the Helmholtz free energy contains 
additional "chemical potential" terms corresponding to the other conserved quantities, Q 
and J, 

F = M -TS -$Q -07 , (5.11) 

where $ = y- and Q — where r + is the Boyer-Lindquist radial co-ordinate at the 
horizon. For a charged hole, the action also contains electromagnetic terms. The surface 
electromagnetic term, eqn. |3i| has the value 4- J d?x\/—h F ab A a rib. However, in order to 
have a regular vector potential, we must gauge transform it to A a — V a (t — $) so that A a 
vanishes on the surface. Hence, the surface action is again given by the gravitational term, 
which has the Euclideanized value Sf UT ^ = —nr + 2 . It is easy to verify using eqn. |5.11| that 



this again leads to a black hole entropy equal to one-fourth of the horizon surface area, and 
an external entropy of zero. 

For non-stationary black holes, the extrinsic curvature also includes a term for the ex- 
pansion of the horizon, K = g + 9. Inserting this into the surface action enables us to 
calculate the instantaneous entropy as matter falls into the membrane in a non-equilibrium 
process. Of course, like the horizon itself, the entropy grows acausally. 

VI. HAMILTONIAN FORMULATION 

The equations of motion can equally well be derived within a Hamiltonian formulation. 
This involves first singling out a global time co-ordinate, t, for the external universe, which 
is then sliced into space-like surfaces, E 4 , of constant t. We can write in the usual way 
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t a = (J^J = aU a -v a , (6.1) 

where U a is the unit normal to E t , U 2 = — 1, and a and — v a are Wheeler's lapse and 
shift, respectively, with v a = the ordinary 3-velocity of a particle with world-line U a . 
For convenience we choose the stretched horizon to be a surface of constant lapse so that 
a, which goes to zero at the true horizon, serves as the stretched horizon regulator. The 
external Hamiltonian for electrodynamics, obtained from the Lagrangian via a Legendre 
transform and written in ordinary three-dimensional vector notation, is 

H out [<p, A, tt] = ^ J d 3 x^g~ Qa(E ■ E + B ■ B) +v- (e x B) - (v • #)) , (6.2) 

where 3 g a b is the 3-metric on S t , <fi = —A a t a is the scalar potential, A a = 3 <7 a 6 Ab is the three- 

— * 

dimensional vector potential, and 7? a = — \rgE a its canonical momentum conjugate. Note 
that E a = F ab U b is the co-moving electric field; E and B above refer to the fields measured 
by a fiducial observer with world-line U a . Finally, the scalar potential is non-dynamical; its 
presence in the Hamiltonian serves to enforce Gauss' law as a constraint. The equations of 
motion are now determined by Hamilton's equations and the constraint: 

SH v 6H u 5H 

-5t =A > lX = -*> H=°- (6 - 3) 
In the bulk these equations are simply Maxwell's equations but, because of the inner bound- 
ary, the usually discarded surface terms that arise during integration by parts now need to 
be canceled. It is easy to show then that the above equations hold only if additional surface 
terms are added to the Hamiltonian: 

H = H out - J d 2 x^ 3s ■ A . (6.4) 

For Maxwell's equations to be satisfied in the bulk, the surface terms are once again the sur- 
face charges and currents necessary to terminate the normal electric and tangential magnetic 
fields at the stretched horizon. Thus, the membrane paradigm is recovered. 

However, it is perhaps more interesting to proceed in a slightly different fashion. In- 
stead of adding new terms, we can use the external Hamiltonian to prove the validity of 

24 



a principle of minimum heat production. Such a principle, which holds under rather gen- 
eral circumstances for stationary dissipative systems, holds for black holes also in slightly 
non-stationary situations. 

Now the time derivative of the external Hamiltonian is not zero, again because of the 
inner boundary. We can use Hamilton's equations to substitute expressions for the time 
derivative of the field and its momentum conjugate. Hamilton's equations are 

A = —aE + vxB — V</> (6.5) 

l = Vx (aB + vx E) , (6.6) 

so that, making repeated use of the vector identity 

V • (V x W) = W • (V x V) - V • (V x W) , (6.7) 

we obtain for the energy loss 
1 f 

- H = -— J d 2 x^ (n ■ (aE\\ xaB^+v- (E ± aE\\ + B ± aB\^) . (6.8) 

So far, we have used only Hamilton's equations. It remains, however, to implement the 
constraint. Hence we may regard —H as a functional of the Lagrange multiplier, (p. We 
therefore have 

That is, the equations of motion follow from minimizing the rate of mass increase of the 
black hole with respect to the scalar potential. This is an exact statement; we now show 
that this reduces to a minimum heat production principle in a quasi-stationary limit. Now 
we note that the first law of black hole thermodynamics allows us to decompose the mass 
change into irreducible and rotational parts: 

dM dQ n dJ , 
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where Qh is the angular velocity at the horizon, and J is the hole's angular momentum. 
Since \v\ — > Qh at the horizon, we see that the second term on the right in eqn. |6.8| 
corresponds to the torquing of the black hole. When this is small, we may approximate 
the mass increase as coming from the first, irreducible term. Hence, in the quasi-stationary 
limit, for a slowly-rotating black hole, the black hole's rate of mass increase is given by the 
dissipation of external energy. Invoking the regularity condition, eqn. |3.12j , then gives 



D[4>] = +^J d2x Vl («4) 2 . JT = > (6- 11 ) 



where atE\\ is given by eqn. |6.5| . This is the principle of minimum heat production: mini- 
mizing the dissipation functional leads to the membrane equation of motion. 

We observe that we could have anticipated this answer. The numerical value of the 
Hamiltonian is the total energy of the system as measured at spatial infinity (assuming an 
asymptotically flat spacetime). The time derivative is then simply the rate, as measured by 
the universal time of distant observers, that energy changes. The rate of decrease of energy 
is the integral of the Poynting flux as measured by local observers, multiplied by two powers 
of a, one power to convert local energy to energy- at- infinity, and one power to convert the 
rate measured by local clocks to the rate measured at infinity. Thus we can immediately 
define a dissipation functional: 

D[(j>) =~J d 2 x^n ■ (E H x B H ) , (6.12) 

where the subscript H denotes that a power of a has been absorbed to renormalize an 
otherwise divergent fiducial quantity. 

In this manner, we can easily write down the dissipation functional for gravity for which 
time-differentiating the Hamiltonian is a much more laborious exercise. The local rate of 
energy transfer is given by the right-hand side of the heat transfer equation, eqn. |3.35| . The 



Hamiltonian for gravity satisfies two constraint equations with the lapse and shift vector 
serving as Lagrange multipliers. Since the membrane picture continues to have a gauge 
freedom associated with time-slicing, the constraint equation associated with the lapse is 
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not implemented. This implies that the dissipation is a functional only of the shift. Hence 
we have 

D[v A \ = J d 2 x^y + 2Wh + a 2 T b a n a U b ) , (6.13) 

where again the two powers of a have been absorbed to render finite the quantities with 
the subscript H. Extremizing D with respect to v leads to the membrane equations of 
motion, enforcing the gauge constraint or, equivalently, obeying the principle of minimum 
heat production. 



VII. CONCLUSION 

We have derived the equations for the membrane paradigm of black holes from an action 
principle directly by demanding that both terms in eqn. |3.1| are stationarized separately. 
This brings advantages of conceptual unity and ease of generalization over the traditional 
approach of manipulating the equations of motion. Specifically, the derivation makes it 
clear why a membrane picture, including dissipative behavior, is possible. A fundamental 
advantage of having an action principle is the guidance it offers for quantization-a property 
we used to fix the constant in the Bekenstein-Hawking formula. 



VIII. APPENDIX 



In this appendix, we shall prove that eqn. |3.21| is zero in the limit that the stretched 
horizon approaches the true horizon. In that limit, an" — > l a . We shall make liberal use 
of Gauss' theorem, the Leibniz rule, and the fact that h ah n^ = K ab rib = 0. In order to use 
Gauss' theorem, we note that since the "acceleration", a c = n d Vd?T, c of the normal vector 
(not to be confused with the fiducial acceleration U d V dU c ) is zero, the 4-covariant divergence 



and the 3-covariant divergence of a vector in the stretched horizon are equal, eqn. [276. 

Now, variations in the metric that are in fact merely gauge transformations can be set 
to zero. Using a vector v a where v a vanishes on the stretched horizon, we can gauge away 
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the variations in the normal direction so that 5g ab — > Shab- Then the left-hand side of eqn. 



3.21 becomes 



J d 3 x^hh bc (V a K<^ c ) - V c (n a 5h ab )) 

= J d 3 xV^h (V a (/i 6c n a 5/i 6c ) - (V a /i fcc ) n a 5/i fec - V c (/i b W/i afe ) + (V c /i bc ) n a 5h, 
= J d 3 xy/^h (v a (h bc n a 5h bc ) + (n c a b + n b a c ) 5h bc - (h hc n a 5h, 



<ab 



-h bc n a 5h ab a c - Kn b n a 5h ab - a b n a 5h ab 



(using h = g — n n c , K ab = +/t£V c n&, and V c w c = w? c + w c a c for w c eH) 
J d 3 xV^h (V« (h bc n a 5h bc ) - Kn b n a 5h, 
(using Gauss' theorem, and a c = 0) 
J d 3 xV^h (v a (h bc %?8h h )j - K (5 (n b n a h 

J d 3 xV^hV a (h hc -l a 5h h )j 
(using h ab n b = 0, and an a — > l a ) 
J d 3 xV^h (h bc -l a 8h, 



'ah 1 - n a h ab Sn b - n b h ab 5n a 



be 



I a 



(8.1) 
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